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Abstract 

Loss of plant diversity influences essential ecosystem processes as aboveground productivity, and can have cascading 
effects on the arthropod communities in adjacent trophic levels. However, few studies have examined how those changes 
in arthropod communities can have additional impacts on ecosystem processes caused by them (e.g. pollination, 
bioturbation, predation, decomposition, herbivory). Therefore, including arthropod effects in predictions of the impact of 
plant diversity loss on such ecosystem processes is an important but little studied piece of information. In a grassland 
biodiversity experiment, we addressed this gap by assessing aboveground decomposer and herbivore communities and 
linking their abundance and diversity to rates of decomposition and herbivory. Path analyses showed that increasing plant 
diversity led to higher abundance and diversity of decomposing arthropods through higher plant biomass. Higher species 
richness of decomposers, in turn, enhanced decomposition. Similarly, species-rich plant communities hosted a higher 
abundance and diversity of herbivores through elevated plant biomass and C:N ratio, leading to higher herbivory rates. 
Integrating trophic interactions into the study of biodiversity effects is required to understand the multiple pathways by 
which biodiversity affects ecosystem functioning. 
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Introduction 

Extensive studies have identified biodiversity loss as one of the 
main drivers of ecosystem change [1], showing that the relation 
between biodiversity and ecosystem processes is mainly positive 
[2,3] . For a long time biodiversity research has focused mainly on 
processes direcdy relating to the trophic level manipulated (e.g. 
aboveground productivity in plant diversity experiments [4,5]; but 
see [6]). More recent studies have shown that plant diversity and 
associated changes in plant species composition also affect 
ecosystem processes governed by higher trophic levels, such as 
decomposition and herbivory [7-10]. 

Decomposition of organic matter is a key process in ecosystems, 
replenishing the pool of plant available soil nutrients and releasing 
photosynthetically fixed carbon back to the atmosphere by the 
activity of several groups of soil organisms. Invertebrate macro- 
fauna, such as isopods, diplopods and earthworms, fragment dead 
plant material and thereby pave the way for further microbial 
decay and mineralization. Consequendy, changes in their density 
and diversity should lead to altered decomposition rates [11,12]. 
They may also affect nutrient cycling and plant productivity of 
grassland ecosystems. Herbivores can decrease plant productivity 
by decreasing plant performance, but contrary to that, they could 
also increase plant productivity by recycling nutrients and 



triggering compensatory growth [13,14]. Arthropod density and 
diversity were discussed as drivers of plant diversity effects on 
decomposition and herbivory, but empirical evidence is scarce 
[8,15]. Studies in experimental grasslands have generally reported 
higher diversity and abundances of herbivores, detritivores, 
pollinators, carnivores and parasitoid arthropods at higher plant 
diversity [16-23], although the strength of the effect differed 
between the groups. 

Plant-animal, plant-ecosystem process, and animal-ecosystem 
process relationships have often been studied separately, but have 
rarely been linked. To understand how plants, arthropods, and 
arthropod-related processes such as decomposition and herbivory 
are linked, there is a strong need to conduct studies considering 
multiple trophic levels and ecosystem processes simultaneously. 
There is a large number of competing hypotheses linking 
consumer-resource dynamics in diverse communities (summarized 
in Fig. 1), applicable to the relationship between arthropod and 
plant communities. Hypotheses emphasizing effects of resources 
on consumers have focused on biomass (More Individuals 
Hypothesis- MIH, Productivity Hypothesis- PH), and stoichiom- 
etry (C:N- Mechanism 3 in Fig. 1), both variables being strongly 
affected by changing plant species richness [24,25] but also on 
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Figure 1. Possible mechanisms driving a positive plant 
diversity-ecosystem process relationship. The consumer-related 
ecosystem processes studied here are herbivory and decomposition. 
Based on existing literature we hypothesize different mechanisms 
driving the positive relationship between plant diversity and herbivory/ 
decomposition, represented here by the different arrows. For details on 
the hypotheses see Introduction. Solid lines indicate expected positive 
correlations, dotted lines expected negative correlations. M1 -3 = Mech- 
anisms 1-3, DM = Dietary Mixing, MIH = More Individuals Hypothesis, 
PH = Productivity Hypothesis, RCH = Resource Concentration Hypothe- 
sis, RHH = Resource Heterogeneity Hypothesis. M1-3 are "unnamed" 
mechanisms, which are reported in the literature. 
doi:10.1371/journal.pone.0106529.g001 

species richness per se (Dietary Mixing- DM, Resource Concen- 
tration Hypothesis- RCH, Resource Heterogeneity Hypothesis- 
RHH). 

MIH and PH hypothesize a positive relationship between plant 
biomass and consumer diversity, but they differ in the underlying 
mechanisms. The MIH [26] applies to communities, which are 
limited in productivity and predicts higher consumer abundances 
in more productive sites as well as consumer diversity as an 
increasing function of total abundance. Following the PH, a higher 
overall resource level in diverse plant communities directly attracts 
more generalist consumer species [27]. The stoichiometry of 
resources is also reported to be a strong predictor of consumer 
abundances as higher food quality (e.g. higher N content of plants) 
increases arthropod fitness (e.g. fecundity) and hence their 
abundance [28]. Following the literature, the way and direction 
how plant species richness affects consumer abundance strongly 
depends on their food specialization. According to the RCH [29] 
we would expect lower abundances of specialized herbivorous 
pests in more diverse habitats due to lower densities of the 
respective host plant. For generalist herbivores the opposite 
pattern may be true (DM): the ability of some generalist herbivores 
(e.g., grasshoppers) to mix their diet to an optimal combination of 
nutrients or to dilute toxins could increase their fitness and hence 
their abundance, indicated by a positive relation between plant 
diversity and arthropod abundance [30]. For consumer species 
richness the RHH [31] predicts that due to a higher number of 
different resources with increasing plant diversity, the number of 
specialist consumer species increases. 

While the theoretical framework of how plant diversity affects 
consumer (arthropod) abundance and diversity is quite elaborated, 



the effects of arthropod communities on herbivory and decompo- 
sition are less clear. Consumer effects on related processes could be 
driven by arthropod abundance (Mechanism 1 in Fig. 1) or species 
richness (Mechanism 2 in Fig. 1). If specialist herbivores dominate, 
then increasing herbivore diversity is likely to increase community- 
wide herbivory rates, while a higher number of herbivores will 
only increase herbivory on those plants that are attacked. By 
contrast, if herbivory is mainly due to generalists, then increasing 
herbivore abundance rather than diversity will increase herbivory. 
For decomposition, functional complementarity among decom- 
posers has been shown in a number of studies [11], yet 
decomposition may well be limited by decomposer abundance. 
Thus, various chains of effects are possible emphasizing the need 
for empirical approaches. All the hypotheses and mechanisms 
(Fig. 1) are not mutually exclusive, however, and many of the 
hypothesized pathways can have interactive and direct or indirect 
effects. The relative importance of these pathways may change 
when considering herbivores as opposed to decomposers because 
some mechanisms like compensatory plant growth would only 
result from herbivory, but not from decomposition. 

Here, we integrate consumers into our models of plant diversity 
effects on ecosystem processes by focusing on three components of 
plant communities (food quantity: aboveground biomass; food 
quality: C :N ratio of plants, and species richness- hereafter: plant 
diversity), two components of consumer communities (arthropod 
abundance, arthropod species richness), and two processes 
(decomposition and herbivory). We test (1) if plant diversity effects 
on arthropods and decomposition/herbivory found in former 
studies can be confirmed, (2) if effects of plant diversity on 
decomposition and herbivory are mediated by altered arthropod 
community composition in terms of overall abundance, and (3) 
which mechanisms link plants to arthropods (plant diversity, 
quantity or quality), and arthropods to decomposition and 
herbivory (abundance or species richness). 

Material and Methods 

Ethic statement 

Arthropod sampling was conducted with the permission of the 
city council of Jena, Germany. 

Experimental field site 

The study was conducted in the framework of The Jena- 
Experiment (Thuringia, Germany, 50°55' N, 11°35' E; 130 m 
above sea level). This grassland biodiversity experiment uses sixty 
species native to regional floodplains. Plant communities were 
sown in 82 plots of 20 x20 m with a gradient of species richness (1, 
2, 4, 8, 16 and 60) and plant functional group richness (1-4 
functional groups: grasses, small herbs, tall herbs, legumes) [32]. In 
2009, the plot size was reduced to 100 m 2 , and two of the 
monocultures were abandoned due to poor cover of the target 
species. The particular species mixtures were randomly selected 
with some constraints [32]. A gradient in abiotic soil properties 
was taken into account by arranging plots in four blocks 
perpendicular to the river Saale [32]. All plots are mown twice a 
year, a management regime common for such meadows. Plots are 
weeded twice a year to maintain the experimental diversity 
gradient. As sown and realized number of target species are highly 
correlated (R 2 = 0.97), we used the number of sown plant species 
for our analyses. The study site was established in May 2002 
allowing for plenty of time for colonization by arthropods and 
making the established plant diversity gradient particularly well 
suited for examining consumer effects on ecosystem functioning 
[20]. 
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Plant community analyses 

In our path analyses we used aboveground plant biomass as 
measure of resource availability for arthropods, because herbivores 
directly feed on plants and decomposers on plant litter material, 
which is reported to be highly correlated with plant biomass [19]. 
In 2010, aboveground plant community biomass was harvested 
during peak standing biomass in late May and August on all plots, 
in the frame of the long-term measurement of aboveground 
productivity. This was done by clipping the vegetation three cm 
above the soil surface in two rectangles of 0.2x0.5 m. The 
harvested biomass was separated into sown species, cleaned from 
weeds, dried at 70°C for 72 h and weighed. As a general measure 
of plant quality we used plant molar C:N ratio, because it is the 
most practicable measure of plant quality across 60 different plant 
species, applicable to many different herbivore and decomposer 
species. Community carbon (C) and nitrogen (N) concentrations of 
the pooled, dried and milled plant material were measured 
separately for May and August 2010 with an elemental analyzer 
(Flash EA 112 Thermo). 

Arthropod sampling 

Arthropods were sampled in 2010 according to methods used in 
previous years with the permission of the city council of Jena. We 
sampled decomposers from May until September 2010, using two 
pitfall traps of 4.5 cm diameter per plot. The use of pit fall traps is 
quite common to sample litter- and surface-dwelling macro- 
arthropods like isopods and diplopods [33]. Traps were emptied 
every two weeks, but closed for two weeks during the mowing 
campaigns in June and September. In total, this resulted in 125 
sampling days over the whole vegetation period. During the 
sampling periods, the field traps were filled with 3% formalin, and 
after emptying the traps, animals were stored in 70% ethanol. The 
resulting activity abundance data will be termed 'abundance' 
hereafter. Our data on decomposers comprise surface active 
decomposer macrofauna only. Less mobile decomposers, such as 
earthworms, may be important agents in surface litter decompo- 
sition [34], but are underrepresented in these samples. 

Herbivorous arthropods were collected by suction sampling in 
June and July using a modified commercial vacuum cleaner 
(Karcher A2500, Karcher GmbH, Winnenden, Germany). Three 
subplots of 0.75 mx0.75 m were randomly chosen within each 
plot, covered quickly by a cage with gauze coat of the same size, 
and sampled until no arthropods were spotted anymore. The 
sampling was carried out between 9 a.m. and 4 p.m. within two 4- 
day sampling periods. Pit fall samples and suction samples were 
sorted to order level and most orders containing decomposers and 
herbivores were further identified to species level. For herbivores 
only Auchenorrhyncha, Coleoptera, and Heteroptera were 
identified to species level, and these orders covered 55% of the 
potential herbivores we sampled. Aphids and Diptera, which made 
together 30.2% of the sampled potential herbivores were not 
identified, because of their patchy distribution within a single plot 
(Aphids) and due to identification difficulties (Diptera). 

Data of the different sampling campaigns were pooled. 

Litter decomposition 

We measured litter decomposition using the standard litterbag 
technique [1 1], We placed one litterbag in each of the 80 plots of 
The Jena-Experiment between 17 th June and 24 th August 2009, as 
successfully done in previous studies [9,10,35]. The fact that we 
have a very high number of replicates within each plant diversity 
level (across plots) compensates for the low replication within plots 
(one litterbag each) since we were more interested in the effects of 
plant diversity than of those of specific plant mixtures. To allow 



access for decomposer macrofauna (e.g. earthworms, isopods, 
diplopods) we constructed litterbags with 4 mm mesh and 
anchored them to the soil surface. We filled each litterbag with 
~3 g of wheat straw (N = 0.4%, C = 45.2%, C:N= 111.5) of 
~3 cm length. At the end of the experiment, the remaining litter 
was cleaned, dried (70°C, 48 h) and weighed to determine mass 
loss. The use of wheat straw acts as a standard method to measure 
potential decomposition rates caused by decomposers and affected 
by microclimatic conditions, without including effects of litter 
quality [9,10,36]. This method focusses on microenvironmental 
effects on decomposition caused by changes in plant communities 
instead of combined effects of microenvironment and litter quality, 
composition or structure. 

Herbivory 

Arthropod herbivory was measured in the plant communities 
during peak biomass in August 2010. A maximum of 30 fully 
developed leaves were sampled randomly for each species from 
each plot where it occurred in the plant biomass samples. The 
minimum number of leaves for a species in a plot was one, and the 
average number was 21 per species and plot. On average 65% 
(SE±2.1%) of the sown species per plot occurred in the biomass 
sampled, therefore the measure can be regarded representative of 
the sown plant communities. The damaged surface area caused by 
herbivores of each leaf (in mm 2 ) was estimated visually by 
comparing the damaged leaf area to circular and square templates 
ranging in size from 1 mm 2 to 500 mm 2 . The leaf area was 
measured with a leaf-area meter (LI-3000C Area Meter, LI-COR 
Biosciences, Lincoln, USA), and the proportional herbivory 
damage (herbivory rates) was calculated as the damaged area 
divided by measured leaf area. The community herbivory rates for 
each plot were calculated by summing the species-specific 
herbivory rates weighted by species-specific leaf biomass in the 
plot. Mammals were excluded from the experimental field site by a 
fence, thus only invertebrates were responsible for the observed 
herbivore damage. Slugs contributed to the measured leaf 
damage, but were not included in herbivore analyses. For details 
in the methods see [8] . 

Statistical analysis 

We analyzed our data using (i) linear models and (ii) path 
analyses. Linear models served as affirmation of the results arising 
from the path analyses, as our sample size (N = 80) was somewhat 
lower than usually recommended for path analyses [37,38]. 

Linear models were used to test for effects of plant community 
properties on arthropods and ecosystem processes. These models 
were fitted using the statistical software package R (R development 
Core Team, http://www.R-project.org. version 2.13.1). Arthro- 
pod abundance, arthropod species richness (decomposer and 
herbivores), decomposition and herbivory were used as response 
variables. For each of the response variables, we fitted a model 
including block (factor), plant diversity (numeric), and functional 
group richness (numeric; number of plant functional groups 
occurring in the respective plot- see experimental design section) as 
explanatory variables. We omitted plant functional group richness 
from our further models, as it had no consistent effect on any of 
the response variables (see Table SI). 

Path analyses were calculated in Amos v.20.0 (IBM SPSS Amos, 
http://www.spss.com/amos/) and were used to infer the hypoth- 
esized causal links between plant diversity and ecosystem 
processes, as ecosystem process responses are likely indirectiy 
mediated by physical, structural or chemical attributes of plants 
and/ or consumer groups studied. We constructed two models (one 
for decomposition and one for herbivory), based on our 
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hypotheses extracted from the published literature (see Introduc- 
tion, and Fig. 1), and did not perform any model simplification. 
Models included arthropod abundance and species richness, 
aboveground plant biomass as a measure of food quantity, and 
plant community C:N ratio as measure of food quality as 
endogenous variables. Because of our low overall sample size 
(N = 80), we did not construct latent variables but assumed direct 
paths among each pair of endogenous variables. As herbivores 
were sampled after the first biomass harvest we used herbivory, 
biomass and C:N data of the data collection in August 2010. For 
decomposition we used mean values of the first and second data 
collection for aboveground biomass and C:N ratio of plants. 

For all analyses (linear models and path analyses) measures of 
plant diversity, decomposer abundance, herbivore abundance and 
decomposer species richness were log transformed, aboveground 
plant biomass was square-root transformed, and herbivory was 
logit transformed [39] to account for non-normality of errors and 
heteroscedasticity, and to linearize the relationships among 
variables. All data used in this publication are deposited in the 
project database (www.the-jena-experiment.de). 

Results 

Effects of plant diversity on the decomposer community 
and decomposition 

Overall, the final dataset contained 1,316 macrofauna decom- 
poser individuals belonging to 12 species (four Isopoda and eight 
Diplopoda species). Isopoda dominated the macrofauna decom- 
poser community, representing on average 92.8% of the total 
community. 

Decomposer abundance and species richness increased signif- 
icantly with increasing plant diversity, but remained unaffected by 
the number of plant functional groups (Table SI, Fig. SI). The 
number of individuals increased from 78.6 per plot (SE±19.2) in 
monocultures to 244 per plot (SE±30.5) in 60-plant-species- 
mixtures. Mean decomposer species richness increased from 3.6 
(SE±0.5) per plot in plant monocultures to 6.3 (SE±0.4) in 60- 
plant-species-mixtures. Decomposition rate increased with in- 
creasing plant species richness from 3 mg mg" 1 d" 1 (SE±0.0002) in 
monocultures up to 5 mg mg -1 d _1 (SE± 0.0004) in 60-plant- 
species-mixtures. The number of plant functional groups, available 
in the respective plant communities (Table SI), did not 
significandy affect decomposition. 

Mechanisms driving altered decomposition rates 

A chi-squared test indicated that our path analysis model cannot 
be rejected as a potential explanation of the observed covariance 
matrix (j? 4 = 4.49, P = 0.344). The model explained 27% of the 
variance in decomposition. The positive effect of plant diversity on 
decomposition (Table S2 and S3, Fig. 2a) was mediated by 
aboveground plant biomass, which increased with increasing plant 
diversity and which positively affected decomposer abundance 
(Fig. 3a) and species richness (marginally). Although decomposer 
abundances and species richness were strongly correlated (Fig. 3b; 
Table SI), only species richness increased litter decomposition 
(Fig. 3c; Table SI). While increasing plant diversity increased 
plant C:N ratio, there was no significant effect of plant C:N ratio 
on decomposers or decomposition. A chi-squared test for the same 
path analysis excluding species richness and abundance of 
decomposers indicated that it need to be rejected as a potential 
explanation of the observed covariance matrix (#%- = 15.60, 
P = 0.016). 



Effects of plant diversity on the herbivore community 
and herbivory rate 

Overall, the final dataset contained 12,829 herbivores, belong- 
ing to 127 species (38 Auchenorrhyncha, 62 Coleoptera, 27 
Heteroptera). On average the herbivore community was domi- 
nated by Coleoptera (48.6%), followed by Auchenorrhyncha 
(43.6%), and Heteroptera (7.8%). Plant diversity, but not plant 
functional group richness, strongly affected the abundance and 
species richness of herbivorous insects and herbivory rate (Table 
SI, Fig. S2). Herbivore abundance increased from 65.2 (SE± 13.7) 
individuals per plot in monocultures to 251 (SE±36.7) individuals 
per plot in 60-plant-species-mixtures. Herbivore communities in 
monocultures contained on average 15.6 (SE±1.0) species per 
plot, whereas in 60-plant-species-mixtures 23.5 (SE±0.8) species 
per plot occurred. Plant diversity also led to a higher herbivory 
rate, nearly doubling from 2% (SE±0.4%) in monocultures to 
3.5% (SE±0.9%) in 60-plant-species-mixtures. 

Mechanisms driving altered herbivory 

A chi-squared test indicated that the path analysis explained 
22% of the variance in herbivory, and cannot be rejected as a 
potential explanation of the observed covariance matrix 
ij 2 4 = 7.50, P = 0.112). The positive relationship between plant 
diversity and herbivory rate (Table S4 and S5, Fig. 2b), was not 
driven by aboveground plant biomass, as there was no significant 
correlation between aboveground biomass and herbivore abun- 
dance and species richness. Herbivore abundance and species 
richness increased significantly with increasing plant diversity (Fig. 
S2). Increasing plant diversity increased plant C:N ratio, which, in 
turn, was positively correlated with herbivore abundance (Fig. 4a). 
Herbivore abundances increased herbivore species richness 
(Fig. 4b), but only higher abundances of herbivorous arthropods 
increased herbivory rate (Fig. 4c). A chi-squared test for the same 
path analysis excluding species richness and abundance of 
herbivores indicated that it need to be rejected as a potential 
explanation of the observed covariance matrix (y 2 6 = 24.70, 
P = 0.001). 

Discussion 

Plant diversity influenced the abundance and diversity of 
decomposers and herbivores, and in turn, these effects propagated 
to decomposition and herbivory rates. Effects of plant diversity on 
process rates could be better explained when including informa- 
tion on the consumer community into the explanatory model for 
both, decomposition and herbivory. Abundance of decomposers 
was positively associated with increased resource (plant) biomass, 
whereas herbivore abundance increased with higher C:N ratio of 
the plant material and plant species richness. However, herbivore 
richness increased directiy with higher plant diversity, whereas 
decomposer richness increased indirectly (via plant biomass and 
decomposer abundance). Accordingly, we provide support for the 
More Individuals Hypothesis (MIH) to drive altered decomposer 
communities, whereas we explain changes in herbivore commu- 
nities by Dietary Mixing (DM) and the Resource Heterogeneity 
Hypothesis (RHH). Decomposer species richness enhanced 
decomposition of litter material (Mechanism 2 in Fig. 1), whereas 
higher herbivory was mainly driven by the higher abundance of 
herbivores (Mechanism 1 in Fig. 1). 

Decomposers and decomposition 

Notably, the discussion of plant community effects (biomass and 
C:N ratio) on the decomposer community is based on to the use of 
standard litter material in our study instead of plot-specific litter 
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Figure 2. Path diagram explaining plant community effects on decomposition and herbivory. Models relate plant community variables 
(diversity, quantity and quality), species richness and abundance of (a) decomposer arthropods to decomposition, and (b) herbivorous insects to 
herbivory. Standardised path coefficients are given on top of the path arrows with significances indicated by *, P<0.05; **, P<0.01; ***, P<0.001. Non- 
significant paths are given in grey. 
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material. Consequently, we here discuss indirect effects between 
plant variables and decomposer communities (plants affect the 
environment for decomposers, which in turn leads to changes in 
the decomposer community) [9,10]. 

Decomposition was positively correlated with plant diversity as 
has been shown in some studies using standard litter [9,40-42] but 
not all [10,34]. Reich and colleagues [3] recently showed that 
plant diversity effects on plant biomass production increased over 
time, and it was suggested that this increase in plant diversity 
effects may be linked to delayed soil feedback effects of slowly 
assembling decomposer communities [43] . Indeed, the density and 
diversity of macro-decomposers and the density of meso-decom- 
posers increased significantly with increasing plant diversity only 
four (aboveground) and six years (belowground) after establish- 
ment of The Jena-Experiment [20,23]. Consequently, the lack of 
plant diversity effects in other studies may be due to the short-term 
character of most previous experiments. Positive effects of plant 
diversity on decomposition processes can possibly be attributed to 
higher productivity of more diverse plant communities and thus 
elevated food availability for decomposers (see MIH in Fig. 1), as 
the amount of available litter material is positively correlated to 
aboveground biomass in grasslands [44]. 

Our path analysis identified aboveground biomass as an 
important driver of changes in the decomposer abundance, which 
supports the MIH. Plant biomass production increased with 
increasing plant diversity and this, in turn, increased decomposer 
species richness via abundance. Resulting increases in decompo- 
sition rates were mediated by increased decomposer species 
richness but not abundance, indicating a more diverse decomposer 
community to be functionally complementary in fragmenting litter 
material (according to M2 in Fig. 1). 

Our results on decomposition are in accordance with a study of 
Heemsbergen et al. [36], who found functional dissimilarity 
among decomposer species being the best predictor for changes in 
leaf litter mass loss and soil respiration. Similar results are reported 
for benthic detritivores in streams, where decomposition rates also 
increase with detritivore richness [45]. In these systems, potential 
mechanisms involve increasing faciliation of resource use among 



species and reduction of inter-specific competiton [46]. Such 
facilitation can arise if the processing of dead organic matter by 
one species enhances the suitability of this resource for another 
species. The positive correlation between decomposer species 
richness and decomposition could also occur due to a sampling 
effect (higher chance of having a functionally dominant decom- 
poser species in a species-rich community than in a species-poor 
one); however, we are not able to disentangle both possible effects 
here. 

Even if we used wheat straw as litter material in our study we 
expected a negative effect of plant C:N ratio on the environment 
for decomposers (lower nutrient availability in litter and soil), 
thereby indirectly affecting their fitness (in our case: abundance). 
In contrast to our hypothesis we detected no significant effect of 
aboveground plant quality (C:N ratio) on decomposer abundance 
and species richness, maybe indicating that the C:N ratio of the 
fresh plant material is a poor predictor of the C:N ratio of the litter 
material. 

In our path analysis a significant direct link between plant 
diversity and decomposition remained, indicating that part of the 
plant diversity effect on decomposition was independent of 
consumer-mediated effects. Potential environmental variables that 
could contribute to this link might include soil moisture of the top 
soil layer, vegetation density, and shade which are important 
predictors of soil dwelling decomposer community composition. 
All three variables are documented to increase with increasing 
plant diversity [10,47]. Further, our analysis was restricted to 
aboveground living decomposers, but we know that belowground 
decomposer groups are strongly affected by plant diversity [19,48]. 
Therefore, we suggest that also belowground decomposition might 
show an imprint of plant species richness via decomposer 
abundance or diversity. As data on soil dwelling decomposer 
species were not available from the the study year and could 
therefore not be included in the analyses this is a natural extension 
for future studies. 
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Figure 3. Pairwise correlations visualizing the significant links 
detected in the path analysis relating plants, decomposers and 
decomposition. We show the relationships between aboveground 
plant biomass and decomposer abundances (a), between decomposer 
abundances and their species richness (b) and decomposer species 
richness and decomposition (c). For statistics, see Table S2. 
doi:10.1371/journal.pone.0106529.g003 

Herbivores and herbivory 

Plant diversity increased both herbivore abundance and 
richness. Plant diversity effects on herbivore richness could not 
be explained by variations in plant biomass or plant C:N, which 
indicates that it may be a direct link between resource diversity 
and consumer diversity. Such relationships are expressly predicted 
for specialist consumers, as more diverse resource assemblages 
enable the presence of a higher number of specialist consumer 
species (e.g. [16,17,21,31]; see RHH in Fig. 1). However, 
generalist consumers can also benefit from more diverse food 
plants, if these provide a nutritionally more balanced or temporally 
less variable food [17,30]. In contrast to findings of McNaughton 
et al. [49] and Borer et al. [22], who found higher aboveground 
plant biomass at higher plant species richness providing the food 
for a larger abundance of herbivores (see MIH in Fig. 1), increases 
in herbivore abundance in this study were not mediated by plant 
biomass. 
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Figure 4. Pairwise correlations visualizing the significant links 
detected in the path analysis relating plants, herbivores and 
herbivory. We show the relationships between plant C:N ratio and 
herbivore abundance (a), herbivore species richness and their 
abundance (b), and between herbivore abundance and herbivory rate 
(c). For statistics, see Table S4. 
doi:10.1371/journal.pone.0106529.g004 

Food plant quality (here C:N ratio) affected the herbivore 
abundance and thereby indirectly herbivory. As previously shown 
[25,50], plant C:N ratios increased with plant species richness. A 
higher food C:N ratio should lead to lower fitness and abundance 
of herbivores [28,51,52] if the plant C:N is above the consumers' 
threshold elemental ratios [53]. Consequendy, community-wide 
herbivory was shown to decrease with poorer autotroph quality 
although individuals compensated for poor food quality by 
ingesting more [54] . If compensatory feeding in a low nutrition 
environment is effective, a neutral relationship between food 
quality and community herbivory might occur [51,55]. Yet, it is 
hard to imagine how compensatory feeding could lead to 
increased herbivore abundance at higher plant C:N as observed 
in our analysis, unless C:N ratio is correlated with other non- 
nutritional aspects of these trophic interactions. Increased 
investment in vertical growth and thus C-rich stem-material in 
the more species-rich plots is a promising candidate for such a 
relationship, which would increase the habitable volume for 
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arthropods and thereby increase abundance. This interpretation is 
supported by the fact that the average height in summer of all 
plants occuring in a communitiy was a significant predictor of 
herbivory levels when modelling community herbivory based on 
plant functional traits [56]. 

Overall, plant diversity increased both herbivore abundance 
and richness, but only herbivore abundance then increased 
herbivory rates (confirming Ml in Fig. 1). The obvious interpre- 
tation of this relationship is that at higher plant diversity a larger 
abundance of herbivores consumes more plant biomass. Conse- 
quently, an increase of herbivory at higher plant diversity [8,57], 
can potentially be explained by an increase in herbivore 
abundance with increasing plant diversity, as documented in this 
study. Integrating data on feeding type or food specialization could 
help for a detailed understanding of the mechanisms responsible 
for the observed positive relation between herbivore abundance 
and herbivory rate. Herbivore species richness was indirecdy 
linked to herbivory rates via its positive correlation to herbivore 
abundance. 

Conclusion 

Our results demonstrate the importance of interactions across 
trophic levels for ecosystem process rates of green and brown food 
webs. Decomposition was driven by arthropod species richness via 
arthropod abundance, and herbivory depended on the number of 
herbivores. Excluding data on consumer levels from our analysis 
(thus measuring only the direcdy plant related links from plant 
diversity and composition to rates of decomposition or herbivory) 
resulted in no explained variance (model rejection, see result 
section) for decomposition and herbivory. In particular, our results 
strongly support the conclusion that integrating trophic interac- 
tions into the study of biodiversity effects is required to understand 
the multiple pathways by which biodiversity affects ecosystem 
functioning [58]. In addition, we could show that the functional 
role of an important part of overall biodiversity - in our case the 
arthropod community - cannot always be deduced a priori, but 
must be unraveled using a mechanistic approach. While functional 
biodiversity research has shown the general importance of 
biodiversity for many ecosystem processes, understanding the 
functional role of biodiversity requires analyzing the mechanisms 
underlying BEF relationships. Given that decomposers and 
herbivores interact with communities of natural enemies exerting 
a top-down control in addition to the bottom-up control by the 
plant community [59] a consequent step forward would be to 
include higher trophic levels (predators, parasites, parasitoids) in 
the analysis of plant diversity effects on arthropod-mediated 
ecosystem processes. 
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